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13.  The  overall  object  of  this  project  was  a  better  theoretical  understanding 
of  the  dynamics  of  energy  rich  molecules  in  both  bound  and  unbound 
states.  Particular  attention  was  to  be  given  to  the  understanding  and 
elucidation  of  selectivity  in  the  role  of  initial  energy  acquisition  and 
its  subsequent  disposal.  Some  effort  was  directed  at  the  understanding  of 
such  selectivity  not  only  for  isolated  molecules  but  also  for 
molecule-surface  collisions.  Two  novel  methodologies  were  developed  and 
broght  to  fruition  during  these  studies.  One  is  the  use  of  algebraic 
Hamiltonians  for  the  description  of  high  overtones  (and  hence  highly 
anharmonic)  energy  rich  states  and  the  determination  of  the  potential 
energy  from  observed  spectra  in  this  way.  Very  extensive  work  has 
demonstrated  the  realistic  feasibility  of  this  approach  for  triatomic 
molecules  (both  linear  and  bent)  and  for  linear  tetratomic  molecules.  A 
complete  summary  of  this  approach  has  just  been  written  up  for  publication 
as  a  book  which  includes  table  of  parameters  of  algebraic  Hamiltonians  for 
a  variety  of  molecules.  The  other  new  theortical  methodology  studies  in 
this  program  is  the  use  of  the  maximum  entropy  information  theory  approach 
for  the  description  of  the  intramolecular  dynamics  of  energy  rich 
molecules  both  below  and  above  the  dissociation  limit.  In  particular, 
this  approach  was  used  to  "invert"  an  observed  Raman  excitation  spectrum, 
that  is,  to  extract  from  the  frequency  spectrum  the  time  evolution  of  the 
system.  Other  techniques,  as  needed,  were  also  used  in  th4  studies  of  the 
time  evolution.  In  particular,  the  use  of  the  algebraic  approach  towards 
the  elucidation  of  the  separation  of  time  scales  and  its  relation  to  the 
restricted  sampling  of  the  available  phase  space  prior  to  dissociation. 
Extensive  runs  of  classical  trajectories  were  also  performed  so  as  to 
demonstrate  the  incomplete  sampling  of  phase  space  even  at  energies  well 
above  dissociation.  The  selectivity  studied  for  surface  processses 
focused  on  the  interplay  between  energetic  and  steric  factors.  It  was 
shown  chat,  as  in  the  gas  phase,  one  can  carefully  di^tnguish  between 
processes  whose  efficiency  either  decreases  dr  increases  with  the 
surface.  The  steric  factors  which  were  discussed  include  the  shape  of  the 
incident  molecule  and  the  'structure  senstivity'  of  the  j process  to  the 
particular  face  of  the  surface. 
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The  overall  object  of  this  project  was  a  better  theoretical  understanding  of  the  dynamics 
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of  energy  rich  molecules  in  both  bound  and  unbound  states.  Particular  attention  was  to 
be  given  to  the  understanding  and  eluddation  of  selectivity  in  the  role  of  initial  energy 
acquisition  and  its  subsequent  disposal.  Some  effort  was  directed  at  the  understanding  of 
such  selectivity  not  only  for  isolated  molecules  but  also  for  molecule-surface  collisions. 
Two  novel  methodologies  were  developed  and  brought  to  fruition  during  these  studies. 
One  is  the  use  of  algebraic  Hamiltonians  for  the  description  of  high  overtones  (and  hence 
highly  anharmonic)  energy  rich  states  and  the  determination  of  the  potential  energy  from 
observed  spectra  in  this  way.  Very  extensive  work  has  demonstrated  the  realistic  feasibil¬ 
ity  of  this  approach  for  triatomic  molecules  (both  linear  and  bent)  and  for  linear 
tetratomic  molecules.  A  complete  summary  of  this  approach  has  just  been  written  up  for 
publication  as  a  book  which  includes  tables  of  parameters  of  algebraic  Hamiltonians  for  a 


variety  of  molecules.  The  other  new  theoretical  methodology  studied  in  this  program  is 
the  use  of  the  maximum  entropy  information  theory  approach  for  the  description  of  the 
intramolecular  dynamics  of  energy  rich  molecules  both  below  and  above  the  dissociation 
limit  In  particular,  this  approach  was  used  to  ’invert’  an  observed  Raman  excitation  spec¬ 
trum,  that  is,  to  extract  from  the  frequency  spectrum  the  time  evolution  of  the  system. 
Other  techniques,  as  needed,  were  also  used  in  the  studies  of  the  time  evolution.  In  partic¬ 
ular,  the  use  of  the  algebraic  approach  towards  the  elucidation  of  the  separation  of  time 
scales  and  its  relation  to  the  restricted  sampling  of  the  available  phase  space  prior  to  dis¬ 
sociation.  Extensive  runs  of  classical  trajectories  were  also  performed  so  as  to  demon¬ 
strate  the  incomplete  sampling  of  phase  space  even  at  energies  well  above  dissociation. 
The  selectivity  studied  for  surface  processes  focused  on  the  interplay  between  energetic 
and  steric  factors.  It  was  shown  that,  as  in  the  gas  phase,  one  can  usefully  distinguish 
between  processes  whose  efficiency  either  decreases  or  increases  with  the  surface.  The 
steric  factors  which  were  discussed  include  the  shape  of  the  incident  molecule  and  the 
’structure  sensitivity’  of  the  process  to  the  particular  face  of  the  surface. 


I.  INTRODUCTION 


The  overall  object  of  this  project  was  a  better  theoretical  understanding  of  the  dynamics 
of  energy  rich  molecules  in  both  bound  and  unbound  states.  Particular  attention  was  to 
be  given  to  the  understanding  and  elucidation  of  selectivity  in  the  role  of  initial  energy 
acquisition  and  its  subsequent  disposal.  Some  effort  was  directed  at  the  understanding  of 
such  selectivity  not  only  for  isolated  molecules  but  also  for  molecule-surface  collisions. 
Two  novel  methodologies  were  developed  and  brought  to  fruition  during  these  studies. 
One  is  the  use  of  algebraic  Hamiltonians  for  the  description  of  high  overtones  (and  hence 
highly  anhamaonic)  energy  rich  states  and  the  determination  of  the  potential  energy  from 
observed  spectra  in  this  way.  Very  extensive  work  has  demonstrated  the  realistic  feasibil- 
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ity  of  this  approach  for  triatomic  molecules  (both  linear  and  bent)  and  for  linear 
tetratomic  molecules.  A  complete  summary  of  this  approach  has  just  been  written  up  for 
publiation  as  a  book  which  includes  tables  of  parameters  of  algebraic  Hamiltonians  for  a 
variety  of  molecules.  The  first  and  last  chapters  of  this  book  are  included  in  Appendix  A. 
The  other  new  theoretical  methodology  studied  in  this  program  is  the  use  of  the  maxi¬ 
mum  entropy  information  theory  approach  for  the  description  of  the  intramolecular 
dynamics  of  energy  rich  molecules  both  below  and  above  the  dissociation  limit.  In  partic¬ 
ular,  this  approach  was  used  to  ’invert’  an  observed  Raman  excitation  spectrum,  that  is,  to 
extract  from  the  frequency  spectrum  the  time  evolution  of  the  system.  A  fully  detailed 
report  on  this  work  is  already  available.  It  is  however  quite  lengthy  so  the  preliminary 
announcement  of  these  results  is  included  as  an  Appendix  B.  Other  techniques,  as 
needed,  were  also  used  in  the  studies  of  the  time  evolution.  In  particular,  the  use  of  the 
algebraic  approach  towards  the  elucidation  of  the  separation  of  time  scales  and  its  rela¬ 
tion  to  the  restricted  sampling  of  the  available  phase  space  prior  to  dissociation.  Exten¬ 
sive  runs  of  classical  trajectories  were  also  performed  so  as  to  demonstrate  the  incom¬ 
plete  sampling  of  phase  space  even  at  energies  well  above  dissociation.  The  selectivity 
studied  for  surface  processes  focused  on  the  interplay  between  energetic  and  steric 


factors.  To  begin  with,  it  was  shown  that,  as  in  the  gas  phase,  one  can  usefully  distinguish' 
between  processes  whose  efficiency  either  decreases  or  increases  with  the  surface.  The 
steric  factors  which  were  discussed  include  the  shape  of  the  incident  molecule  and  the 
’structure  sensitivity’  of  the  process  to  the  particular  face  of  the  surface. 

We  turn  now  to  a  discussion  of  specific  results.  The  reference  numbers  are  to  the  list 
of  publications.  Section  6  of  this  report. 

2.  Algebraic  Hamiltonian 

At  the  conclusion  of  this  proposal  it  is  possible  to  state  that  the  algebraic  approach 
has  been  fully  validated  for  triatomics.  This  includes  not  only  all  known  vibration- 
vibration  coupling  but  also  vibration-rotation  coupling.  The  transcription  from  an  alge¬ 
braic  to  a  geometrical  picture  is  fully  at  hand  [4,  Ch.  7  of  31]  and  the  connection  with 
intramolecular  dynamics  is  being  explored  [26].  Our  plans  for  the  future  in  that  direction 
are  outlined  in  Appendix  A. 

3.  Intramolecular  Dynamics 

Three  reviews  [1,18,25]  provide  a  summary  of  our  work  up  to  a  year  ago.  These  are 
quite  detailed  and  cover  the  entire  gamut  of  our  activities  so  it  would  be  superfluous  to 
repear  it  and  we  include  the  latest  of  those  as  Appendix  C.  Instead,  we  draw  attention  to 
the  most  recent  work  (for  a  short  report,  see  Appendix  B)  where  the  maximum  entropy 
formalism  has  been  effectively  employed  to  represent  spectra.  Not  discussed  therein  are 
implications  towards  a  separation  of  time  scales.  This  is  shown  in  figure  8.1  of  Appendix 
A.  The  connection  to  our  previous  AFC3R  supported  work  on  intensity  fluctuations  has 
been  firmly  established  and  is  in  preparation  for  publication.  The  different  aspects  of 
intramolecular  dynamics  as  discussed  in  the  original  proposal  are  all  extensively  repre¬ 
sented  in  the  list  of  publications. 
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4.  Sticking  and  Dissociation  in  Molecule-Surface  Collision 

Following  the  experimental  demonstration,  by  the  late  R.B.  Bernstein,  of  a  large 
steric  effect  in  the  sacking  of  molecules  to  surfaces  and  taking  cognizance  of  the  AFOSR 
interest  in  selectivity  of  surface  processes,  we  attempted  a  theoretical  discussion  of  this 
effect  [6].  The  theory  clearly  predicted  a  decrease  of  the  effect  with  increasing  collision 
energy  with  the  surface.  We  therefore  checked  the  theory  for  dissociation  processes 

t 

whose  probability  is  known  to  rise  steeply  past  the  energetic  threshold  [13].  When  this 
could  be  explained  we  examined  the  sensitivity  to  the  structure  of  the  surface  [15,23]. 
Work  is  continuing  in  this  direction,  supported  by  a  PRF  grant  to  Dr.  FaiTelly,  who  is  by 
now  in  Utah. 

A 

5.  Cooperation  with  Other  AFOSR  Projects 

Two  papers  are  a  result  of  discussions  with  other  AFOSR  supported  scientists,  dis¬ 
cussions  which  took  place  at  the  annual  contractors  meetings.  Paper  [7]  is  a  discussion  of 
a  nossible  experimental  test  of  the  extent  of  geometrical  randomization  prior  to  dissocia¬ 
tion.  Papers  [21]  and  [30]  discuss  the  extent  of  energy  randomization  following  selective 
excitation.  The  groups  of  Professors  Thompson  and  Raff  have  since  done  more  work  in 
that  direction. 
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Appendix  A: 


ALGEBRAIC  THEORY  OF  MOLECULES 


This  appendix  consists  of  the  introduction  and  the  concluding  chapter  from  a  book 
of  this  title  (by  F.  Iachello  and  R.D.  Levine)  which  will  be  published  in  1993  by  Oxford 
University  Press. 

The  entire  ms.,  which  is  over  500  pages,  is  available  upon  request. 


i 

I 


INTRODUCTION 


Molecular  ::  iscopy  is  undergoing  an  essential  change.  Seemingly  the  change  is 

only  quantitative;  better  initial  state  preparation,  improved  light  sources  and  specially 
designed  pumping  schemes  and  more  sensitive  detection  techniques  are  providing  ever 
improved  resolution  and  a  wider  range  of  accessible  final  states.  Figure  I. 

Figure  0.1.  Stimulated  emission  pumping  (SEP,  Hamilton  et  ah,  1986;  Northrup 
and  Sears,  1992)  is  a  new  experimental  technique  for  accessing  higher  lying  vibra¬ 
tional  levels  of  molecules  in  their  ground  electronic  states.  Shown  is  the  SEP  vibra¬ 
tional  spectrum  of  S02  where  a  pair  of  dips  represent  one  vibrational  level. 
(Adapted  from  Yamanouchi,  Takeuchi  and  Tsuchiya,  1990).  The  stick  spectrum  at 
the  bottom  represents  the  position  of  the  vibrational  levels  given  by  Equation  0.1 
with  the  constants  given  in  Table  0.1.  The  bright  levels  are  represnted  by  longer 
sticks. 

A  closer  examination  suggests  also  a  qualitative  change.  New  ideas  and  not  only  better 
results  are  forthcoming.  One  example  of  the  changing  attitudes  is  the  increasing  concern 
with  time  evolution.  The  time  -  energy  uncertainty  relation  and  the  pursuit  of  higher  res¬ 
olution  means  that  traditional  spectroscopy  implicitly  is  equivalent  to  the  study  of  the  sta¬ 
tionary  states  determined  by  the  long  time  limit  of  the  intramolecular  dynamics.  The 
recent  increasing  concern  with  the  role  of  anharmonicities  and  of  resonance  couplings, 
made  unavoidable  by  the  study  of  higher  lying  rovibrational  states  and  the  experimental 
reality  of  avoiding  inhomogenous  broadening  (e.g.,  using  supersonic  expansion.  Quack, 
1990)  make  the  entire  time  domain  of  direct  interest  to  spectroscopists  (Bitto  and  Huber, 
1992).  The  very  complementarity  with  the  studies  in  the  frequency  domain  (broad  homo¬ 
geneous  spectral  features  h  early  time  dynamics  and  vice  versa)  makes  lower  resolution 
spectra  of  interest.  On  the  other  hand,  the  traditional  concerns  of  spectroscopy  (Herzberg, 
1945,1950;  Barrow,  1962;  King,  1964;  Hollenberg,  1970;  Herzberg,  1971;  Bunker,  1979; 
Steinfeld,  1985)  remain  very  much  with  us.  Figure  0.2. 
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Figurt-  0.2.  Direct  ovenone  spectroscopy  of  C2H2  using  Fourier  transform  spec¬ 
troscopy.  Here,  at  high  resolution,  the  entire  band  of  rotational  transitions,  which 
accompany  a  given  vibrational  transition,  can  be  resolved.  Here  the  band,  in  the  vis¬ 
ible  range,  corresponding  to  the  direct  excitation  of  v  =  5  of  the  v3  stretch  mode  is 
shown.  (Adapted  from  Herman  et  al.,  1991.  See  also  Scherer,  Lehmann  and  Klem¬ 
perer,  1983  and  Figure  8.4). 

Better  determination  of  overtone  and  combination  bands  of  familiar  molecules  and  the 
spectroscopic  characterization  of  new  species  (radicals  (Shida,  1991;  Bemath,  1990),  ions 
(Miller  and  Bondybey,  1983;  Leach,  1980)  and  van  der  Waals  molecules  in  particular, 
Nesbitt,  1988;  Saykally,  1989;  Hutson,  1990;  Heaven,  1992)  continue  to  receive  wide 
attention.  Figure  0.3. 

Figure  0.3.  The  structure  of  the  van  der  Waals  molecule  QHu  •  Ar  as  determined  by 
very  high  resolution  spectroscopy.  (Adapted  from  Weber,  van  Bargen,  Riedle  and 
Neusser,  1990).  The  potential  along  the  QH*  -  Ar  stretch  motion  is  shown  in  Figure 
1.4. 

What  is  needed  for  modem  spectroscopy  is  a  formalism  able  to  discuss  both  level 
structure  beyond  the  harmonic  limit  and  the  corresponding  dynamics.  A  Hamiltonian  is 
thus  unavoidable  for  it  is  the  generator  of  time  evolution.  Yet  there  needs  to  be  a  practical 
method  for  the  determination  of  the  eigenvalues  of  this  Hamiltonian.  As  in  the  traditional 
Dunham-like  expansion,  it  will  be  useful  if  the  spectra  can  be  well-approximated  by  a 
small  number  of  constants.  An  example  of  such  an  expansion  is  the  fit  of  the  vibration 
spectrum  of  S02  in  its  electronic  ground  state,  Figure  0.1,  by  the  expansion 

G(V| ,  v2,  v3)  =  £  <y,(v,  +  1/2)  +  £  Xjj(v,  +  1/2 )(v,  +  1/2) 

« 

+  X  y,jk(v,  +  1/2 )(V;  +  l/2)(v*  +  1/2)  , 

iS/S* 


(0.1) 
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Here  <7(vt,  v2,  v3)  is  the  level  energy  in  wave  number  units  (as  far  as  possible  we  follow 
the  notation  of  Herzberg,  1950)  and  the  constants  in  Equation  0.1  are  given  in  Table  0.1. 

Table  0.1.  Vibrational  constants'^  (cm-1)  for  the  ground  electronic  state  of  S02 


o>\ 

1167.84(15) 

Jin 

-0.0061(10) 

oh 

522.21(19) 

>222 

-0.0014(18) 

oh 

1382.18(50) 

>333 

-0.031(22) 

-3.655(23) 

>112 

-0.0001(17) 

*22 

-0.374(36) 

>113 

-0.1574(41) 

*33 

-5.36(20) 

>122 

-0.0063(19) 

*12 

-3.129(29) 

>123 

-0.0509(53) 

*13 

-14.277(71) 

>133 

0.255(11) 

*23 

-4.122(8) 

>223 

0.0214(50) 

>233 

-0.008(12) 

fl)A  number  in  parentheses  represents  a  standard  error  obtained  by  the  least-squares 
analysis.  The  vibrational  constants,  <u„  xtj,  and  yijt  are  the  expansion  coefficients  of 
Eq.  (0.1).  Adapted  from  Yamagouchi,  Takeuchi  and  Tsuchiya,  1990. 

As  usual  the  v’s  are  the  vibrational  quantum  numbers  of  S02  and  rather  high  (above  10) 
values  can  be  reached  using  the  SEP  technique.  Eq.  (0.1)  provides  a  fit  to  the  observed 
levels  to  within  an  error  below  10  cm-1  which  is  almost  the  experimental  accuracy.  We 
need  however  to  be  able  to  relate  the  parameters  in  this  expansion  directly  to  a  Hamilto¬ 
nian.  The  familiar  way  of  doing  this  proceeds  in  two  steps.  First,  the  electronic  problem 
is  solved,  in  the  Bom-Oppenheimer  approximation,  leading  to  the  potential  for  the 
motion  of  the  nuclei.  Then  the  Schrddinger  equation  for  the  eigenvalues  of  this  potential 
is  solved.  Since  for  any  but  diatomic  molecules,  the  potential  is  a  function  of  many  coor¬ 
dinates,  neither  the  first  nor  the  second  steps  are  simple  to  implement.  For  a  number  of 
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test  cases  this  procedure  has  been  carried  out  and  for  diatomic  molecules  of  lower  rows 
atoms  it  can  challenge  experiment  in  its  precision.  For  larger  molecules  it  is  still  not 
practical  to  compute  the  required  potential  with  sufficient  accuracy.  It  is  therefore  often 
approximated  using  convenient  functional  forms.  Not  too  far  from  a  deep  equilibrium 
point,  the  potential  can  be  expanded  in  the  displacements  coordinates  relative  to  the  equi¬ 
librium  configuration.  Such  a  ’force  field’  representation  is  quite  convenient  but  is  of 
limited  validity  for  higher  lying  states  due  to  the  slow  convergence  of  such  a  power-series 
expansion.  More  flexible  functional  forms  that  can  describe  the  asymptotic  dissociation 
plateaux  (Murrell  et  al.,  1984)  require  many  parameters.  The  purely  numerical  solution 
of  the  SchiOdinger  differential  equation  for  the  eigenvalues  of  such  a  potential  makes  the 
optimization  of  the  parameters  in  the  potential,  via  a  fit  to  the  observed  spectrum,  a  large 
scale  numerical  problem  complicated,  as  all  such  non-linear  problems  are,  by  non  unique¬ 
ness  and  by  local  minima. 

In  this  book  we  present  an  alternative  approach.  Our  discussion  in  this  introductory 
volume  will  put  particular  emphasis  on  the  traditional  concerns,  namely  determing  the 
levels  and  intensities  of  the  corresponding  transitions.  The  approach  we  present  retains, 
at  least  in  part,  the  simplicity  of  a  Dunham-like  approach  in  that,  at  least  approxiately,  it 
provides  the  energy  as  an  analytic  function  of  the  quantum  numbers  as  in  Eq.  (0.1).  If 
this  approximation  is  not  sufficient,  the  method  provides  corrections  derived  in  a  system¬ 
atic  fashion.  On  the  other  hand,  the  method  starts  with  a  Hamiltonian  so  that  one  obtains 
not  only  eigenvalues  but  also  eigenfunctions.  1  It  is  for  this  reason  that  it  can  provide 
intensities  and  other  matrix  elements. 

The  Hamiltonian  used  in  our  approach  is  an  algebraic  one  and  so  are  all  the  opera¬ 
tions  in  the  method.  This  unlike  the  more  familial,  differential  operators  of  wave  mechan¬ 
ics.  The  technical  advantage  of  an  algebraic  approach  is  the  comparative  ease  of  alge¬ 
braic  operations.  Equally  important  however  is  the  result,  obtained  by  comparison  with 
experiment,  that  there  are  generic  forms  of  'Jgebraic  Hamiltonians  and  that  entire  classes 
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of  molecules  can  be  described  by  a  common  Hamiltonian  where  only  the  (typically,  lin¬ 
ear)  parameters  are  different  for  the  different  molecules. 

The  algebraic  (or  matrix)  formulation  of  quantum  mechanics,  footnote  1,  is  less 
familiar  than  the  differential  (or  wave)  formulation.  This  is  a  disadvantage  and  one  pur¬ 
pose  of  the  present  volume  is  to  show,  by  explicit  examples,  the  benefits  of  the  algebraic 
approach.  The  interested  reader  will  have  to  judge  if  the  benefits  are  sufficient  to  over¬ 
come  the  potential  barrier  to  the  understanding  of  a  new  approach.  We  intend  to  demon¬ 
strate  that  the  algebraic  formulation  is  indeed  a  viable  alternative. 

The  algebras  one  uses  are  Lie  algebras.  These  algebras  were  introduced  at  the  end 
of  the  19th  Century  by  Lie  but  it  is  only  in  the  much  more  recent  past,  footnote  2,  that 

A 

they  are  being  used  in  physics.  In  the  approach  that  we  follow,  emphasis  is  put  on  start¬ 
ing  with  the  Hamiltonian  and  the  main  technical  tool  is  the  algebra  rather  than  on  the  cor¬ 
responding  group.  In  the  appendices  to  this  book,  we  give  an  account  of  the  important 
properties  of  Lie  algebras  that  are  of  interest  for  the  applications  in  the  study  of 
molecules,  footnote  3. 

Even  if  one  restricts  attention  to  vibrations  and  rotations  of  molecules,  there  are  a 
variety  of  Lie  algebras  that  one  can  use.  In  some  applications,  the  algebras  associated 
with  the  harmonic  oscillator  are  used.  We  mention  these  briefly  in  Chapter  1.  We  prefer 
however  already  in  zeroth  order  to  use  algebras  associated  with  anharmonic  oscillators. 
Since  an  understanding  of  the  algebraic  methods  requires  a  comparison  with  more  tradi¬ 
tional  methods,  we  present  in  several  parts  of  the  book  a  direct  comparison  with  both  the 
Dunham  expansion  and  the  solution  of  the  Schitfdinger  equation. 

In  this  book  we  shall  write  the  Hamiltonian  as  an  (algebraic)  operator  using  the 
appropriate  Lie  algebra.  We  intend  to  illustrate  by  many  applications  what  we  mean  by 
this  cryptic  statement.  It  is  important  to  already  emphasize  that  one  way  to  represent  such 
a  Hamiltonian  is  as  a  matrix.  In  this  connection  we  draw  attention  to  one  important  area 
of  spectroscopy,  that  of  electronically  excited  states  of  larger  molecules,  footnote  4, 
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which  is  traditionally  discussed  in  terms  of  matrix  Hamiltonians  the  simplest  of  which  is 
the,  so  called,  picket  fence  model  (Bixon  and  Jortner  1968).  A  central  issue  in  this  area  of 
spectroscopy  is  the  time  evolution  of  an  initially  prepared  non  stationary  states.  We  defer 
a  detailed  discussed  of  such  topics  to  a  sequel  volume  which  deals  with  the  algebraic 
approach  to  dynamics. 

In  this  book  we  deal  mainly  with  stationary  states,  their  energies  and  matrix  ele¬ 
ments.  Unless  otherwise  stated,  we  use  the  wave  number  (cm-1)  as  a  measure  of  the 
energy.  The  conversion  factors  with  other  units  are  shown  in  Table  0.2. 

Table  0.2.  Energy  conversion  factors"* . 


Unit 

cm  1 

eV 

Hz(sec-1) 

J/mole 

°K 

1  cm-1 

1 

1.2398(-4) 

2.997924(10) 

1.1962(1) 

1.4388 

leV 

8.06573(3) 

1 

2.41804(14) 

9.6487(4) 

1.16049(4) 

1  Hz(sec-1) 

3.3356(-ll) 

4.1356(-15) 

1 

3.9903(-lU) 

4.7993(-ll) 

1  cal/mole 

0.34976 

4,3363(-5) 

1.04854(10) 

4.1840 

0.50325 

1°K 

0.69503 

8.6170(-5) 

2.0836(10) 

8.1340 

1 

a)Thc  numbers  in  parentheses  represent  powers  of  ten. 

lA  =  10-8cm  =  0. 1  nm  =  10-4  fim  =  108/(cm-1).  Speed  of  light  in  vacuum 
c  =  2.99792458  x  108  msec-1.  The  meter  is  now  defined  as  the  distance  traveled  by 
light  in  vacuum  in  1/c  sec,  so  that  c  is  a  defined  constant  and  the  second  (as  mea¬ 
sured  by  atomic  clocks)  is  the  fundamental  unit.  The  cal  is  still  an  often  used  energy 
unit.  For  a  complete  account,  see  Mills  (1988). 

The  present  volume  deals  exclusively  with  rotation-vibration  spectra.  Electronic 
excitations  can  also  be  described  algebraically,  but  this  description  is  still  at  a  too  prelimi¬ 
nary  stage  for  inclusion  in  a  book  format. 


Footnotes 


1.  The  matrix  formulation  of  quantum  mechanics  was  introduced  as  early  as  1925  by 
Heisenberg,  1925;  Bom  and  Jordan,  1925;  Dirac,  1925;  Bom,  Heisenberg  and  Jordan, 
1926. 

2.  It  was  only  in  the  1930’s  that  Lie  algebras  were  being  used  in  physics  (Weyl,  1931; 
Wigner,  1931,  1937;  van  der  Waerden,  1931;  Yamanouchi,  1937;  Racah,  1942,  1949). 
Most  of  the  early  applications  dealt  with  the  algebra  of  rotations  (Wigner-Racah  algebra). 
In  this  volume  we  aim  to  present  a  unified  approach  which  starts  from  the  algebra  as  the 
key  tool  for  the  construction  of  spectra  and  intensities.  This  approach  originated  in  the 
1960’s  (Dothan,  Gell’Mann  and  Be’eman,  1965;  Barot  and  Btihm,  1965)  and  had  major 
applications  in  nuclear  physics  (Iachello  and  Arima,  1974;  Arima  and  Iacheilo,  1975.  For 
an  update,  see  the  books  by  Iachello  and  Arima,  1987  and  Iachello  and  van  Isecker,  1991 
and  the  reviews  in  Bohm,  Ne’eman  and  Barutr,  1988).  In  our  discussion  we  shall  empha¬ 
size  the  role  of  the  algebra  rather  than  that  of  the  corresponding  group.  Specifically,  we 
shall  seek  such  a  description  that  even  in  zeroth  order  the  spectrum  is  anharmonic.  We 
shall  also  try  to  forge  a  clear  link  with  the  geometrical  point  of  view.  There  are  many 
other  important  applications  of  Lie  groups  and  algebras.  In  particular,  we  do  not  discuss 
time  dependent  aspects  (Aihassid  and  Levine,  1977;  Wulfman,  1979;  Levine  and  Wulf- 
man,  1979;  Levine,  1985)  nor  other  applications  (e.g..  Coherent  States  (Klauder  and 
Skagerstam,  1985;  Pcrelemov,  1986;  Zhang  et  ah,  1990)  or  configuration  interaction  in 
electronic  structure  calculations  (Judd,  1968;  Paldus,  1978;  Hinze,  1981;  Pauncz  and 
Matsen,  1986;  Adams  et  al.,  1987;  Duch,  1988))  where  the  group  structure  is  of  central 
importance. 

3.  Lie  groups  and  Lie  algebras  are  discussed  in  many  textbooks  (Hamermesh,  1962;  Lip- 
kin,  1965;  Gilmore,  1974;  Wyboume,  1974;  Barut  and  Rqczka,  1986).  We  follow  closely 
the  notation  of  Wyboume,  1974.  There  are  also  a  number  of  mathematical  texts  (Miller, 

X 

1968;  Talman,  1968;  Vilenkin,  1968;  Miller,  1977,  Olver,  1986). 
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4.  Hamiltonians  expressed  in  matrix  forms  have  been  extensively  employed  in  the  theory 
of  radiationless  transitions  of  electronically  excited  sates  of  larger  molecules  (Bixon  and 
Jortner,  1968;  Schlag  et  ah,  1971;  Freed,  1972;  Nitzan  et  al.,  1973;  Avouris  et  al.,  1978; 
Jortner  and  Levine,  1981;  Felker  and  Zewail,  1988;  Seel  and  Domcke,  1991). 
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CHAPTER  8 


PROLOGUE  TO  THE  FUTURE 

At  this  point,  we  hope  to  have  demonstrated  that  the  algebraic  approach  provides  a 
viable  method  for  the  quantitative  description  of  molecular  vibrotational  spectra.  Chap¬ 
ters  5  (triatomic  molecules,  both  linear  and  bent)  and  5  (linear  tetratomic  molecules)  and 
Appendix  C  provide  extensive  documentation  for  the  quantitative  applications,  while  Ch. 
6  shows  that  larger  molecules  can  also  be  treated.  Throughout,  but  most  particularly  in 
Ch.  7,  we  have  sought  to  forge  link  with  the  more  familiar  geometrical  approach.  It  is 
precisely  our  requirement  that  even  in  zeroth  order  the  Hamiltonian  we  start  with 
describes  an  anharmonic  motion,  which  makes  this  link  not  trivial.  The  advantage  of  our 

A 

requirement  in  providing,  even  in  zeroth  order,  high  overtone  spectra  which  are  typically 
more  accurate  than  10  cm-1,  should  not  be  overlooked.  Yet  much  remains  to  be  done.  In 
this  chapter  we  look  to  the  future:  Where  and  why  do  we  think  that  the  algebraic 
approach  will  prove  particularly  advantageous.  Of  course,  what  we  really  hope  for  is  to 
be  surprised  by  unexpected  new  developments  and  applications.  Here,  however,  is  where 
we  are  certain  that  some  of  the  future  progress  will  be  made,  with  special  reference  to  the 
spectroscopy  of  higher  energy  states  of  molecules. 

One  area  of  spectroscopy  where  the  Hamiltonian  in  matrix  form  is  the  route  of 
choice  is  that  of  large  polyatomics,  particularly  so  when  in  an  electronically  excited  state, 
footnote  0.3.  Such  states  are  isoenergetic  with  very  high  vibrational  overtones  of  the 
ground  electronic  states  so  that  a  fully  geometrical  approach  is  impractical.  Even  at  lower 
energies,  the  exceedingly  high  density  of  vibrational  states  strongly  favours  an  alternative 
approach  and  the  use  of  model  Hamilton  matrices  is  not  uncommon.  Such  model  matrices 
are  introduced  so  as  to  account  for  the  regularities  that  often  survive  in  the  observed  spec¬ 
trum.  One  such  striking  feature  is  often  referred  to  as  a  ’clump’  (Hamilton  et  al.  1986). 
Consider  a  pure  vibrational  progression  of  states  as  can  be  observed  in  stimulated  emis¬ 
sion  spectroscopy,  SEP,  Figure  0.1.  A  clump  is  a  spectral  feature  that  appears  as  a  broad 


line  at  a  given  spectral  resolution  but  which  is  revealed  to  be  a  set  of  overlapping  lines 
when  a  higher  resolution  is  used.  To  qualify  as  a  ’clump’,  the  width  needs  to  be  inherent 
rather  than  instrumental.  By  inspection  of  Figure  0.1  one  can  imagine  that  the  spectrum 
will  exhibit  clumps  at  a  lower  resolution.  It  can  even  be  that  there  are  clumps  within  a 
clump.  Figure  8.1. 

Figure  8.1.  A  spectral  clump  and  its  two  tiers  of  fine  structure.  (The  sharp  feature  at 
-8770  cm-1  in  panel  b  is  the  envelope  of  the  set  of  lines  shown  in  panel  c).  When 
states  are  coupled  in  a  sequential  manner  there  can  often  be  a  separation  of  time 
scales  in  the  temporal  evolution  of  a  non  stationary  state  or  a  corresponding  separa¬ 
tion  of  frequency  scales  in  the  spectrum.  (Remac  '  and  Levine,  1992).  Algebraic 
Hamiltonians  provide  a  convenient  framework  for  the  discussion  of  sequential  cou¬ 
pling. 

For  small  polyatomics  (e.g.,  N02,  C2FQ  one  can  think  of  the  coupling  of  different 
electronic  states  from  a  geometric  point  of  view  (Englman,  1972  ;  Kb'ppel  et  al.,  1984; 
Whetten  et  al.,  1985).  Even  then  a  transformation  to  a  matrix  Hamiltonian  is  typically 
very  useful  and  is  often  implemented  particularly  when  only  a  very  limited  number  (two, 
three)  electronic  states  are  coupled.  By  using  a  Hamiltonian  expressed  in  terms  of  genera¬ 
tors  for  the  electronic  degrees  of  freedom  one  can  describe  a  multistate  electronic  spec¬ 
trum  (Frank  et  al.,  1986,  1989,  1991;  Lemus  et  al.,  1992;  Lemus  and  Frank,  1991,  1992). 
Moreover,  one  can  couple  the  electronic  motion  to  anharmonic  motion  of  the  nuclei.  In 
the  future  one  can  confidently  expect  to  see  many  more  such  applications. 

An  algebraic  approach  to  electronic  degrees  of  freedom  will  be  particularly  advanta¬ 
geous  when  the  density  of  electronic  states  is  high  and  the  coupling  to  the  iovibrational 
motion  is  strong.  Such  is  the  case  for  very  high  Rydberg  states  which  lie  just  below  the 
ionization  continuum.  For  larger  molecules  there  is  a  very  facile  energy  exchange 
between  the  electronic  and  nuclear  manifold  of  states.  There  is  currently  a  revival  of 
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interest  in  such  problems  (Schlag  and  Levine,  1992).  For  smaller  molecules  (NO,  H20) 
the  relevant  couplings  are  currendy  being  delineated  by  spectroscopic  techniques  (e.g., 
Bryant  et  al.,  1992;  Gilbert  and  Child,  1991). 

Next,  we  turn  to  dynamics  proper.  Spectroscopy  itself  is  intimately  related  to 
intramolecular  dynamics,  that  is  to  the  time  evolution  of  a  non  stationary  state  of  the 
Hamiltonian.  Detailed  applications  of  algebraic  methods  to  such  problems  have  so  far 
been  mostly  limited  to  problems  where  the  Hamiltonian  is  a  linear  expression  in  terms  of 
the  generators  of  the  algebra.  The  reason  why  this  is  the  case  can  be  understood  even 
without  an  appeal  to  the  group  concept:  Let  X,  at  time  t  =  0  be  a  generator  whose  time 
evolution  X(t)  is  of  interest.  Then,  in  the  Heisenberg  picture:  iHdX(t)/dt  =  [X(t),  H).  It 
follows  from  the  closure  property  of  Lie  algebras,  Equation  2.1,  that  if  X(t)  and  H  are  in 
the  algebra  so  is  the  time  rate  of  change  of  X(t)  and  therefore  so  is  X{t  +  dt).  Hence,  if  X 
at  t  -  0  is  an  operator  in  the  algebra  and  the  Hamiltonian  is  any  linear  combination  of 
generators,  X(t)  remains  in  the  algebra.  A  key  open  problem,  footnote  1,  is  how  to  pro¬ 
ceed  in  the  more  realistic  case  when  the  Hamiltonian  is  bilinear  rather  than  linear  in  the 
generators  of  the  algebra. 

The  second  set  of  problems  in  dynamics  are  those  of  scattering  theory  where  the 
Hamiltonian  is  of  the  form  H  =  HQ  +  V  and  the  interaction  V  vanishes  when  the  colliding 
particles  are  far  apart.  It  is  usually  assumed  that  the  H0  part  is  already  solved  and  that  the 
interesting  or  the  hard  part  is  to  account  for  the  role  of  V.  For  realistic  systems  which  arel 
anharmonic  even  the  role  of  H0  can  be  quite  significant  An  example  which  received 
much  recent  attention  is  the  reaction  of  vibrationally  excited  HOD  with  H  atoms  (Sinha  et  ' 
al.,  1991),  Figure  8.2. 

Figure  8.2.  Schematic  illustration  of  the  preferential  reactivity  of  HOD  with  an  H 

atom  when  the  stretch  excitation  is  localized  in  either  one  of  the  two  bonds. 

The  large  difference  in  the  OH  and  OD  vibrational  frequencies  means  that  the  stretch 
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overtones  of  HOD  are  primarily  local  in  character,  cf.  Section  4.21.  It  follows  that  one 
can  excite  HOD  to  overtones  localized  preferentially  on  either  one  of  the  two  bonds  and 
that  an  approaching  H  atom  will  abstract,  preferentially  either  an  H  or  a  D  atom.  Figure 
8.2.  Even  in  the  lowest  stretch  states  this  selectivity  is  maintained  (Bronikowski  et  al., 
1991). 

As  in  scattering  theory  in  general,  one  can  treat  the  role  of  V  in  either  a  time  inde¬ 
pendent  or  a  time  dependent  point  of  view.  The  latter  is  simpler  if  the  perturbation  V  is 
either  explicitly  time  dependent  or  can  be  approximated  as  such,  say  by  replacing  the 
approach  motion  during  the  collision  by  a  classical  path.  Algebraic  methods  have  been 
particularly  useful  in  that  context,  footnote  2,  where  an  important  aspect  is  the  description 

A 

of  a  realistic  level  structure  for  H0.  Figure  8.3  is  a  very  recent  application  to  electron- 
molecule  scattering. 

Figure  8.3.  The  vibrational  elastic  and  inelastic  differential  cross  sections  for  elec¬ 
tron  scattering  off  LiF  at  £=5.44  eV.  (Alhassid  and  Shao,  1992).  Solid  lines:  with 
an  improved  dipole  interaction  (which  breaks  the  0(4 )  symmetry).  Long  dashed 
lines:  the  calculations  by  Bijker  and  Amado  (1986).  The  short  dashed  lines  are  the 
Bom  approximation. 

In  time  dependent  collision  theory,  the  asymptotic  states  ’after’  or  ’before’  the  colli¬ 
sion  are  defined  by  the  limits  of  t  — »  ±oo.  In  time  independent  theory  one  can  either  use  a 
Lippmann-Schwinger  like  approach  (Alhassid  and  Levine,  1985)  or  define  scattering 
through  a  coordinate  representation  (Alhassid  et  al.,  1983, 1984,  1986;  Ojha,  1986;  Wu  et 
al.,  1987,  1989;  Frank  and  Wolf,  1984;  Frank  et  al.,  1986,  1985).  The  latter  has  been 
more  extensively  studied,  since  for  scattering  by  a  potential  the  geometrical  connection  is 
well  understood,  cf.  Chapter  7.  Applications  to  systems  with  internal  structure  (e.g., 
Wehrhahn,  1990,  1991;  Wehrhahn  and  Levine,  1991)  are  still  in  their  early  stages. 

In  its  initial  stage  of  development,  the  algebraic  approach  has  sought  to  show  why 
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and  how  it  provides  a  framework  for  die  understanding  of  large  amplitude,  anharmonic 
motion.  Much  of  our  concern  in  this  volume  has  been  with  this  stage  including,  of  course, 
the  even  simpler  problem  of  uncoupled  motions.  In  its  current  stage  the  approach  will 
seek  to  keep  pace  with  the  new  developments  of  modem  spectroscopy  including  the 
interest  in  floppy  molecules.  The  better  understanding  of  the  overtones  of  polyatomics 
including  the  higher  energy  regimes  where  the  non  rigidity  is  sufficient  to  allow  for  iso¬ 
merization  will,  no  doubt,  continue  to  be  of  central  concern.  This  will  include  not  only 
vibradon-rotation  coupling  but  also  rovibrarion-electronic  coupling.  In  terms  of  the  con¬ 
nection  to  the  geometrical  picture  this  means,  of  course,  that  one  is  interested  in  the  shape 
of  the  potential  further  and  further  away  from  its  global  minimum  and,  in  particular,  in 

A 

the  possibility  of  local  minima.  Moreover,  at  increasing  coupling,  the  very  concept  of  the 
potential  becomes  less  useful.  The  algebraic  approach  offers  a  special  advantage  under 
such  circumstances. 

At  higher  light  intensities,  it  is  no  longer  possible  to  treat  the  coupling  to  the  field  in 
lowest  order  in  perturbation  theory  (George  et  al.t  1977)  and  algebraic  techniques  again 
become  particularly  useful.  The  emphasis  of  a  time  domain  interpretation  of  spectroscopy 
is  very  likely  to  further  increase.  This  is  not  only  because  of  spectroscopic  experiments  in 
the  time  domain  (Khundkar  and  Zewail,  1990;  Pollard  and  Mathies,  1992)  but  also 
because  of  the  ’time  window’  which  is  revealed  by  experiments  in  the  frequency  domain 
(Imre  et  al.,  1984),  a  time  window  which  covers  the  very  range  of  interest  to  chemists. 
The  emphasis  on  the  time  domain  has  also  revived  the  interest  in  coherence  effects  in 
optical  spectroscopy  (Steinfeld,  1985).  Then  there  are  larger  molecules  up  to  and  includ¬ 
ing  molecules  of  biological  interest,  clusters,  etc.  One  has  both  experimental  and  theoreti¬ 
cal  reasons  to  expect  that  increasing  the  size  does  not  necessarily  imply  an  exponential 
increase  in  the  complexity  and  unwieldiness  of  the  spectrum.  Figure  8.4. 

Figure  8.4.  The  v  =  6  CH  overtone  spectrum  as  determined  by  photoacoustic 

absorption  (Hall,  1984)  for  increasing  size  alkvnes.  Each  panel  contains  the 
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molecular  formula  and  the  density  of  vibrational  states  per  cm"1.  The  band  origin 

does  shift  some  with  size  but  is  roughly  at  18450  cm"1  for  a’’  the  molecules  shown. 

(See  also  Kerstel  et  al.f  1991). 

Finally,  one  expects  spectroscopy  to  have  an  increasing  impact  not  only  on 
intramolecular  but  also  on  intermolecular  dynamics.  This  includes  not  only  a  better 
understanding  of  the  reactants,  as  in  Figure  8.2,  but  also  the  application  of  spectroscopy 
to  the  very  collision  process  (Brooks,  1988;  Neumark,  1992).  At  lower  light  intensities, 
spectroscopy  acts  in  its  usual  role  as  a  probe  except  that  it  is  probing  very  transitory 
events.  At  higher  intensities  one  can  guide  and  not  only  monitor  the  evolution  of  the  sys¬ 
tem  (Jortner  et  a!.,  1991;  Rabitz  and  Shi  1991;  Rice,  1992;  Brumer  and  Shapiro,  1992). 
The  discussion  of  control  of  quantum  systems  is  one  more  area  to  which  algebraic  tech¬ 
niques  can  and  will  usefully  contribute. 

i 

i 

! 
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1.  The  problem  of  time  evolution  for  a  Hamiltonian  bilinear  in  the  generators  (Levine, 
1982)  has  been  extensively  discussed.  The  proposed  solutions  .nclude  the  use  of  varia¬ 
tional  principles  (Tishby  and  Levine,  1984),  mean  field  self  consistent  methods  (Meyer  et 
al.,  1988),  time  dependent  constants  of  the  motion  (Levine,  1982)  and  numerous  others 
which  we  hope  to  discuss  in  detail  in  a  sequel  to  this  volume. 

2.  Algebraic  methods  have  been  effectively  employed  in  time  dependent  collision  theory 

(Alhassid  and  Levine,  1977,  1978).  In  the  sudden  approximation  (Levine  and  Wulfman, 

1979)  and  the  use  of  an  H0  which  is  bilinear  in  the  generators  presents  no  real  problems. 

This  has  been  very  effectively  employed  in  electron-molecule  scattering  (Bijker  et  al., 

/* 

1986,  1988,  1990;  Ginocchio  et  al.,  1986;  Mengoni  and  Shirai,  1988,  1991;  Alhassid  et 
al.,  1993),  where  due  to  the  high  velocity  of  the  incident  electron  one  can  neglect  the 
change  in  its  path  due  to  energy  transfer  to  the  me  lecule. 
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Appendix  B: 


THE  DETERMINATION  OF  TIME  CROSSCORRELATION  FUNCTIONS 
BY  INVERSION  OF  RAMAN  EXCITATION  PROFILES 


This  appendix  is  the  entire  text  of  the  short  paper,  with  the  above  title,  by  F.  Remade, 
R.D.  Levine  and  J.L.  Kinsey. 

A  much  more  detailed  version  is  also  ready  for  publication. 


The  Determination  of  Time  Crosscorrelation  Functions 
by  Inversion  of  Raman  Excitation  Profiles 


F.  Remade*  and  R.D.  Levine 

The  Fritz  Haber  Research  Center  for  Molecular  Dynamics 
The  Hebrew  University  of  Jerusalem 
Jerusalem  91904,  Israel 

and 

J.L.  Kinsey 

Rice  Quantum  Institute  and 
Department  of  Chemistry 
Rice  University,  Houston,  TX  77251,  U.S.A. 


The  short  time  dynamics  of  a  nonstationary  state  is  here  determined  by  a  direct 
inversion  of  an  experimentally  measurable  Raman  excitation  profile.  The  procedure  is 
based  on  the  analytic  properties  of  the  Raman  scattering  amplitude  and  is  best  imple¬ 
mented  within  the  maximum  entropy  formalism.  A  computational  example  for  resonance 
Raman  excitation  profiles  in  the  B  state  of  iodobenzene  is  presented. 


PACS  numbers:  33.10.Cs,  33.20.Fb,  02.70.-k1,  33.70.-w 


Chargfi  de  Recherches,  FNRS,  Departement  de  Chimie,  University  de  Liege,  B-4000  Liege  1, 
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Rabinovitch  and  coworkers.  Chemical  clocking5-11-12  in  particular  provides  striking  evidence 
for  exceedingly  fast  delocalization.  Some  of  the  reported  rates  are  just  an  order  of  magni¬ 
tude  slower  than  typical  vibrational  frequencies. 

Our  question  is  whether  the  loss  of  the  selectivity  of  the  initial  excitation,  in  highly 
excited  molecules,  is  indeed  so  rapid.  We  shall  bring  forth  arguments  in  support  of  a  much 
more  gradual  and  sequential  loss.  In  particular  we  shall  argue  that  the  rapid  sub-ps  initial 
stage  will  have  the  same  temporal  character  in  a  purely  harmonic  molecule.  At  the  end  of 
this  first  stage  the  molecule  has  by  far  not  sampled  the  majority  of  its  energetically  accessi¬ 
ble  phase  space.  Hence  we  shall  conclude  that  the  fast  initial  dephasing  is  definitely  not 
sufficient  to  establish  the  representative  coverage  of  phase  space  which  is  assumed  in  the 
RRKM  approach. 

The  reason  for  reexamining  this  question  are  the  recent  high  resolution  spectroscopic 
experiments  of  higher  overtones  of  polyatomic  molecules12-20..  Adopting  a  time  dependent 
point  of  view21-22  such  high  resolution  experiments  provide  insight  on  the  temporal  evolu¬ 
tion  of  the  initially  localized  excitation22-30.  These,  and  the  corresponding  studies  in  the  fre¬ 
quency  domain31-36  furnish  what,  from  the  standpoint  of  chemical  kinetics,  is  complemen¬ 
tary  the  short  time  information. 

The  information  provided  by  both  collisional  clocking  and  the  inversion  of  spectro¬ 
scopic  data  is  essentially  indirect  information.  Ultimately,  real  time  monitoring37-40  will 
allow  a  direct  probe  of  the  dissociation  process.  For  the  time  being,  overtone  spectroscopy 
provides  useful  insights  on  the  first  dozen  ps  or  so  following  the  t  =  0  initiation  of  a  non  sta¬ 
tionary  molecular  state.  It  would  be  of  considerable  interest  to  extend  the  time  domain  into 
the  sub-ns  range.  This  however  is  unlikely  to  be  realized  in  the  immediate  future.  Such  time 
intervals  are  long  on  the  spectroscopic  scale  as  they  require  a  frequency  resolution  which  is 
typically  unrealistic  at  the  high  excitation  energies  of  interest.  Computations  are  equally 
better  at  shorter  times,  for  both  numerical  reasons  and  for  the  need  to  know  finer  and  finer 
details  of  the  Hamiltonian  in  order  to  obtain  a  higher  resolution  spectrum  or,  equivalently, 
an  accurate  longer  time  propagation. 

THE  EXPLORATION  OF  PHASE  SPACE 

Of  the  many  possible  spectroscopic  techniques  we  here  select  for  discussion  the  two 
that  most  readily  correspond  to  the  information  of  interest.  The  normalized  absorption  spec¬ 
trum  S(a}) 22  and  the  normalized  resonant  Raman  excitation  spectrum  tr(at)41,42.  (Note  that 
there  are  other  Raman  cross  sections,  most  notably  the  Raman  emission  cross  section, 
which  is  not  here  discussed).  The  two  types  of  spectra  are  related  to  two  time  correlation 
functions.  The  absorption  spectrum  is  the  Fourier  transform  of  the  time  autocorrelation 
function  C(r)  of  an  initial  non-stationary  state  0(0) 

C(t)  =  J  exp  {-icot)S(a})do} 


-  J  - 


C(r)  =  <*(0MO>  • 

The  spectral  intensity  pH,  at  the  energy  of  the  n’th  eigenstate 

S(Ol)-XPn^~^n)  . 

H 

is  related  to  <fi(0)  by 

Pn  =  I  <  >  I2  • 

The  time  autocorrelation  function  determines  the  'survival  probability’,  IC(r)l2  of  the  non 
stationary  state,  IC(r)l2  =  I  <  $>(/)l^(0)  >  I2. 

It  is  equally  of  interest  to  determine  not  only  where  the  initial  state  is  not,  but  where  it 
is.  This  information  is  provided  by  the  time  crosscorrelation  function  <  $>/(0)l$>,(r)  >  which 
is  related  to  the  resonant  Raman  excitation  cross  section 

oo 

CTfi(a))-=  I J  <  <p f\<pi(t)  >  exp(frar)drl2  . 
o 

/• 

Note  however  that  unlike  S(co)  and  C(r)  which  are  a  Fourier  transform  pair,  the  time  cross 
correlation  function  determines  the  Raman  excitation  cross  section  but  the  inversion  of  this 
relation  is  less  immediate.  The  direct  route  is  thus  to  compute  the  time  crosscorrelation 
function  from  the  dynamics  and  to  obtain  the  excitation  profile  as  the  one-sided  Fourier 
transform.  Recently  we  have  however  developed  a  method  of  analysis  for  extracting  the 
time  crosscorrelation  function  from  the  Raman  excitation  profile.  C(t )  on  the  other  hand 
can  either  be  computed  or  determined  directly  from  the  experimental  data,  figure  1. 


Figure  1.  The  survival  probability  IC(r)i2  determined  by  a  Fourier  transform  of  the 
excitation  spectrum  of  acetone  in  the  visible  and  near  UV.  High  resolution  experimen¬ 
tal  data  of  H.  Zuckermann,  Y.  Haas,  M.  Drnbbels,  J.  Heinze,  L.  Meerts,  J.  Reuss,  and 
J.  van  Bladel. 


/  .  .  . 
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The  short  time  evolution  of  C(r)  has  been  extensively  discussed22-30.  In  particular, 
considerable  attention  was  given  to  the  notion  of  the  correlation  hole  •  •  .  Here,  how¬ 

ever,  we  draw  attention  to  a  much  more  prominent  feature,  namely  the  very  rapid  drop  of 
IC(r)l2  from  an  initial  value  of  unity  to  an  (average)  value  smaller  by  cne  or  more  orders  of 
magnitude,  figure  1.  The  point  is  that  this  drop  occurs  on  a  time  scale  of  vibrational  periods 
and  does  not  provide  evidence  for  intramolecular  vibrational  energy  redistribution.  In  other 
words,  consider  two  Hamiltonians  H0  (whose  spectrum  is  regular)  and  H  -HQ  +  V  where 
V  includes  the  diverse  anharmonic  couplings.  While  V  will  cause  krai  level  repulsions, 
splittings,  etc.,  the  envelope  of  the  spectrum  of  H0  and  H  will  be  essenti  Uy  identical.  Since 
the  gross  structure  in  5(tu)  is  the  same,  the  short  time  dependence  of  C(t)  will  be  the  same. 

That  the  initial  time  evolution  does  not  discern  the  fine  spectral  details  (which  indicate 
the  presence  or  absence  of  chaos)  is  not  a  uniquely  quantal  phenomena.  It  can  just  as  well 
be  understood  on  purely  classical  considerations.  For  historical  reasons  we  phrase  the  dis¬ 
cussion  in  terms  of  Slater’s  model  of  unimolecular  dissociation44.  Recall  that  in  the  Slater 
model  the  molecule  is  perfectly  harmonic.  The  extension  of  the  bond  to  be  broken  is  written 
as  a  linear  superposition  of  normal  modes.  Each  mode  with  its  own  phase,  which  varies 
with  time  in  accord  with  its  harmonic  frequency.  As  the  phases  vary,  so  does  the  extension 
of  the  bond  to  be  broken.  Dissociation  is  said  to  take  place  when,  due  to  a  p.iase  matching, 
the  bond  is  extended  to  a  critical  value.  Now  let  us  run  this  mechanical  system  backwards  in 
time:  An  initially  extended  bond  rapidly  returns  to  a  more  average  value  due  to  the  dephas¬ 
ing  of  the  normal  modes,  each  beating  with  its  own  frequency.  The  analogy  wv.ii  the  optical 
excitation  is  established  by  thinking  of  the  initial  bond  extension  as  resulting  from  a  vertical 
Condon  transition.  The  initial  drop  of  lC(r)l2  corresponds  to  the  Slater-like  dephasing  pro¬ 
cess  of  the  non- interacting  modes. 

In  a  purely  harmonic  Slater-type  model,  the  initial  dephasing  leading  to  the  contraction 
of  the  extended  bond  will  eventually  lead  to  a  rephasing  and  an  elongation  of  the  bond. 
There  is  however  an  important  proviso  since  real  molecules  are  anharmonic.  IC(r)l2  does 
not  recur  back  to  its  initial  value  of  unity  but  to  a  lower  value.  It  is  this,  much  more  gradual, 
relaxation  of  IC(r)l2  which  reflects  the  sampling  of  phase  space.  The  same  is  true  for  a  quan¬ 
tal  time  evolution. 

ANALYSIS  OF  SEQUENTIAL  EVOLUTION  IN  PHASE  SPACE 

The  molecular  eigenstates  provide  one  suitable  basis  for  discussing  the  time  evolution 
of  C(r).  In  that  basis  the  time  history  is  entirely  expressed  in  terms  of  dephasing  since  the 
eigenstates  themselves  are  stationary 

C(t)  =Hpn  exp {-iEntth)  . 

It  is  the  distribution  of  the  intensities,  pm,  and  the  spacings  En  -  En,  as  reflected  in  the 


spectral  autocorrelation  function  G(oj),  the  so  called  'power  spectrum’*3 

IC(f)l2  =  J  dct)cxp(-ia)t)G(co ) 

G(co)  =  £  Z  P„PmS{con  -  com  -  CO) 

n  m 

that  determines  the  different  time  regimes  in  IC(r)l2  if  such  regimes  exist. 

G(cj)  is  the  weighted  distribution  (by  pnp„)  of  spacings  in  the  spectrum.  Note  that 
G{co)  describes  not  only  the  distribution  of  spacings45  but  also  the  fluctuation  of  intensi¬ 
ties.46  Our  understanding  of  molecular  spectra  is  sufficient  for  us  to  recognize  that  G(<y) 
may  well  contain  several  peaks.  The  one  at  highest  spacings  (or  shortest  times)  corresponds 
to  spacings  of  different  vibrational  states.  Next  are  contributions  from  the  dominant  anhar- 
monicitics  (the  largest  resonant  couplings),  etc.  One  sees  very  clearly  the  existence  of  sev¬ 
eral  scales  in  G(o)  for  spectra  taken  at  different  resolutions  (sometimes  referred  to  as 
’clumps’,  see,  e.g.  Refs.  15,18,47).  At  each  increasing  scale  there  is  a  new  set  of  lines, 
which  merge  into  a  clump  when  the  resolution  is  reduced,  figure  2. 
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Figure  2.  A  (computed)  spectrum  at  three  levels  of  resolution,  plotted  vs.  the  frequency 
measured  in  units  of  the  mean  spacings  between  eigenstates.  A  group  of  resolved  lines 
at  the  highest  resolution  (bottom  panel)  appears  as  a  single  feature  (a  ’clump’)  at  a 
lower  resolution.  A  group  of  clumps  appears  as  a  single  clump  at  the  lowest  resolution 
(uppermost  panel). 
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The  nested  structure  of  the  observed  spectrum  suggests  that  it  may  be  useful  to  con¬ 
sider  an  alternative  basis  set,  corresponding  to  a  zeroth  order  Hamiltonian  H q.  The  initially 
excited  optically  ’bright’  states  correspond  to  a  small  subset  of  states  of  H0.  Elsewhere48  we 

have  discussed  the  nested  structure  of  the  perturbation,  V,V  —  Vl+V2-i - ,  which  gives 

rise  to  a  sequence  of  Hamiltonians  Hit  H:  =  H0  +  Vj,  H2  =  Hi  +  V2, ....  such  that  it  is  the 
eigenstates  of  //,•  which  span  the  phase  space  sampled  by  the  /’th  stage.  In  terms  of  the 
eigenstates  of  H0,  each  stage  corresponds  to  a  larger  subset  of  states  which  are  coupled  to 
the  initially  prepared  state.  One  can  therefore  define  projection  operators  M„  which  project 
onto  the  subspace  of  states  of  H0  relevant  at  the  f  th  stage,  such  that 

A/,-  =  +Mi  =  2  . 

Here,  Mr  is  the  projection  on  the  new  set  of  states  brought  into  coupling  at  the  r’th  stage 
and  M0  =  M0. 

The  experimental  observation  of  a  clump  structure  in  the  spectrum  means  that  one  can 
compute  a  survival  probability  for  a  given  clump.  Therefore  one  can  determine  from  the 
observed  spectrum  not  only  IC(r)i2, 

IC(r)l2  =  Tr(M0/y(r))  , 

where  pit)  =  l^(r)  ><  pit) I,  but  also 

IC,(r)l2  =  Tr(Mip(t)) 

=  IW/Ip(;))  . 

« 

In  the  maximum  entropy  approach, 27,48  we  determine  pit)  as  the  most  spread  out  den¬ 
sity  subject  to  the  constraints  imposed  by  the  spectrum,  i.e.  the  values  of  the  survival  prob¬ 
abilities  IC,(r)l2,  t  =  0, 1, . . ..  Taking  into  consideration  that  7>(p(r))  =  1  and  A?w  M(  =  0,  one 
obtains  for  p(t) 

p{t)  =  exp(2  pr(t)Mr)/Z 

r=0 

where  Z  insures  the  normalization  of  pit) 

Z  =  7>(exp(£  prit)Mr)) 

r=0 

R 

and  R  is  the  highest  level  of  nesting  evident  in  the  spectrum,  A//?+I  =I-%Mr  and 

r* 0 

MR+iit)  =  0.  Two  very  suggestive  forms  of  pit),  obtained  using  the  orthogonality  of  the  M' s 

.  *+1 

pit)  =  Z"  2  exp [prit))Mr 

r  r*0 


arc 


1 


/ 
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showing  that  exp (pr{t))Tr[M  r]/Z  is  the  fractional  occupancy  of  the  r’th  subspace  ar.d  the 
nested  form 

p(t)  =  Z~l{MR+l  +  exp(A«(f)) (Mg  +  exp(A/?_1(r))/'iV/s_1  +  ••• 

+  exp(2  i  (t)){M i  +  exp(AQ(f))A/o/  •••}  . 

Here  the  Ar(r)’s  are  the  Lagrange  multipliers  conjugate  to  the  survival  probabilities  !Cr(r)l2 
and  are  related  to  the  pri.t)'s  by 

t-R  -  Mr 
Ar-i  =  Mr- i  “  Mr 

=  Mo  -Mi 

A  plot  of  the  entropy  of  p(t )  for  a  system  manifesting  a  separation  of  time  scales  is  shown 
in  figure  3.  (The  corresponding  frequency  spectrum  is  shown  in  figure  2).  A  detailed  dis¬ 
cussion  of  sequential  time  evolution  is  available  in  Ref.  48.  Earlier  closely  relevant  work 
includes  Refs.  26,  27  and  30. 


o' 

C- 

c 
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Figure  3.  The  entropy  (in  units  of  R )  vs.  time  (in  units  of  ft/ mean  spacings,  logarithmic 
scale)  computed  for  the  density  in  phase  space  of  maximal  entropy.  The  survival  prob¬ 
abilities  vs.  time  are  those  determined  by  the  spectrum  shown  in  figure  2.  Note  the 
clear  separation  into  three  time  regimes  (the  logarithmic  abscissa  scale  makes  it  possi¬ 
ble  to  show  the  different  regimes  on  the  same  plot).  IV  is  the  number  of  eigenstates  so 
that  for  a  strictly  uniform  covering  of  phase  space  the  long  time  value  of  the  entropy 
should  be  R  In  //.  As  can  be  seen  the  actual  asymptotic  value  is  somewhat  lower.  This 
is  due  to  the  fluctuation  of  the  spectral  intensities.'44 

In  a  wider  context,  the  separation  of  time  scales,  as  reflected  in  the  sequential  evolu¬ 
tion  stages  in  phase  space,  is  a  particular  case  of  impulsive  coupling.  Impulsive  coupling, 
manifested  in  a  variety  of  intramolecular  dynamics  systems49  e.g.,  van  der  Waals  clusters. 
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energy  rich  molecules,  reactions  in  liquids50  and  on  surfaces,51  is  the  rather  sudden  .switch 
from  a  motion  confined  to  a  smaller  region  of  phase  space  to  a  kcs  confined  mode.  It  results 
from  terms  in  the  Hamiltonian  which  are  ’non-resonant’  on  a  shorter  time  scale  and  thus  not 
included  for  shorter  time  dynamics. 

DYNAMICS 

An  alternative  to  extracting  the  time  correlation  functions  from  the  observed  spectra  is 
to  compute  them.  This  is  particularly  attractive  because  it  is  the  earlier  time  dynamics  which 
determine  the  more  gross  features  of  the  spectrum.  The  considerable  interest  in  time- 
dependent  methods  as  reflected  in  these  proceedings  is  no  doubt  related  to  this  inherent 
advantage.  Our  own  interest  is  however  in  the  dynamics  as  such  with  special  relevance  to 
unbound  systems.  The  possibility  of  mode  selective  chemistry,32  as  mentioned  in  the  intro¬ 
duction,  is  an  important  motivation.  Here,  we  review  two  recent  studies,  one  in  the  gas 
phase53  and  one  in  solution,54  which  were  based  on  the  use  of  classical  dynamics:  An 
ensemble  of  initial  conditions  which  represents  the  initially  prepared  norv stationary  state,  is 
generated.  Trajectories  are  then  propagated  and  the  time  correlation  function  is  computed. 

Four  molecules  at  energies  both  below  and  above  their  dissociation  threshold  were 
studied  in  the  gas  phase53:  CH3ONO,  C2H4,  SiF4  and  the  large  (by  our  standards)  trinicro  s- 
triazine  (RDX).  Initial  excitations  included  both  mode  selective  and  random  deposition  of 
the  energy.  The  computations  were  run  for  10  ps,  which  is  short  for  thermal  unimolecular 
reactions  where  the  excess  energy  is  low.  We  have  however  examined  all  four  molecules 
also  at  an  energy  where  their  lifetimes  for  dissociation  are  comparable  or  shorter  than  10  ps. 

The  central  conclusion  of  this  work,  parts  of  which  are  still  in  progress,  is  that  on  the 
timescale  of  10  ps  mode  mixing  is  very  incomplete.  This  is  reflected  both  in  the  behavior  of 
the  time  correlation  functions  and  in  the  resulting  spectra.  Figure  4  is  a  typical  outcome  for 
SiF4  at  an  excess  energy  of  75.8  kcal/mole.  The  time  autocorrelation  of  the  SiF  stretch 
modes  even  at  longer  times  is  still  distinct  from  that  of  the  FSiF  bends.  Many  more  results 
are  available  in  the  full  report.53 

For  all  molecules  there  is  a  very  clear  initial  drop  of  IC(r)l2.  As  expected  from  the  ear¬ 
lier  discussion  of  the  Slater  model,  this  drop  is  most  extreme  for  the  many  mode  RDX 
molecule.  However,  this  drop  (quite  moderate  for  SiF4,  figure  4),  does  not  represent  any 
loss  of  selectivity  of  the  initial  excitation.  Taking  figure  4a  as  an  example,  the  initial  drop  of 
IC(r)l  means  that  a  single  SiF  stretch  is  not  a  good  zero  order  mode.  Rather,  it  should  be 
better  thought  of  as  a  linear  superposition  of  the  stretch  normal  modes.  These  rapidly 
dephase  with  respect  to  one  another,  but  the  energy  remains  in  the  stretches.  Plucking  an 
SiF  stretch  shows  no  spectral  activity  in  the  bend  region.  The  same  incomplete  energy 
scrambling  is  seen  even  well  above  the  dissociation  threshold,  figure  5.  (Note  that  the 
energy  of  SiF4  is  high  enough  for  dissociation  to  take  place  in  less  than  10  ps.). 
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The  onset  of  mode  mixing  occurs  later  and  precludes  the  full  recurrence  of  IC(r)l.  It  is 
important  to  note  that  in  these  polyatomics,  the  recurrence  time  can  be  quite  long.  A  simple 
estimate  is  as  follows.  Say  the  initially  prepared  state  is  a  stretched  bond.  Recurrence  means 
that  this  bond  is  again  stretched  to  the  same  extent.  The  mean  time  between  two  such  events 
is  just  what  is  computed  in  the  Slater  theory:  In  that  model,  the  dissociation  lifetime  of  a 
molecule  equated  to  this  mean  recurrence  time. 

SiF4  is  a  prototype  of  molecules  which  have  been  excited  by  infrared  multiphoton 
absorption.  Such  experiments  have  typically  tended  to  support  extensive  energy  scrambling 
in  the  pumped  molecule.1,55  Our  work  suggests  that  this  conclusion  may  well  reflect  the 
long  duration  of  the  typical  C02  lasers  used  in  these  experiments.  (At  least  in  the  ns  scale  if 
not  longer).  It  should  therefore  be  of  interest  to  return  to  such  experiments  but  with  pulse 
lengths  in  the  ps  range. 

The  computations  in  the  liquid54  studied  the  time  evolution  of  HOC1  pumped  to  the 
v  =  6  overtone  of  OH.  This  provides  sufficient  energy  to  dissociate  the  weak  OC1  bond.  Yet 
there  are  good  reasons  to  expect  the  very  non-resonant  energy  flow  from  tjie  OH  to  the  OC1 
stretch  to  be  quite  hindered56  and  to  proceed  primarily  via  the  bend.  The  novel  feature,  very 
evident  in  the  simulations,54  is  the  significant  enhancement  of  this  intramolecular  energy 
transfer  by  collisions  with  the  surrounding  Ar  liquid.  This  was  particularly  evident  when 
the  spectrum  S(co)  is  computed  not  over  the  entire  time  history  (as  it  should  be)  but  over  a 
narrower  time  window.  This  results,  of  course,  in  a  broad  spectrum,  devoid  of  fine  structure. 
However,  the  three  vibrational  modes  of  HOC1  are  so  different  in  their  frequency  that  for  a 
time  window  of,  say,  1  ps,  one  can  readily  discern  the  three  (two  stretches  and  the  bend) 
spectral  regions.  Such  studies  clearly  demonstrate  that  the  intramolecular  vibrational  energy 
redistribution  is  slow  enough  that  the  collisions  with  the  atoms  of  the  IL  fid  can  exercise  a 
beneficial  effect.  From  an  overall  kinetic  point  of  view,  the  density  of  the  liquid  has  addi¬ 
tional  roles.54  It  can  cage  the  dissociation  products.  It  can  also  remove  energy  by  collisions 
with  the  activated  HOC1  molecule.  Even  so,  at  low  densities,  the  primary  role  of  the  sur¬ 
rounding  Ar  atoms  is  to  enhance  the  intramolecular  energy  transfer. 


CONCLUDING  REMARKS 

As  can  be  seen  from  the  results  presented  here,  the  selectivity  of  the  initial  excitatii 
may  well  persist  on  time  scales  comparable  to  the  lifetime  of  an  unbound  molecule.  This 
means  that  the  system  has  not  enough  time  to  sample  all  the  energetically  accessible  phase 
space  before  dissociation  occurs. 

In  highly  excited  polyatomic  molecules,  the  recurrence  time  tr  =  Hp(E)  (p(E)  is  the 
density  of  states),  for  which  the  sampling  of  phase  space  is  by  definition  completed,  is  very 
long.  From  the  standard  expression57  of  the  RRKM  rate  constant  k(E) 


k(E)  = 


N*(E-Eq) 


v>J  * 


»/p(£) 


I 
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(N*(E  -  E0 )  is  the  number  of  states  of  the  transition  state  and  E0  the  threshold  energy),  one 
can  see  that  the  fraction  of  molecules  exp (~N*(E  -  E0))  that  have  not  dissociated  by  r,  is 
very  small  since  N*(E  -  E0 )  is  typically  very  large.  It  is  therefore  not  complete  sampling  of 
phase  space  which  needs  to  be  assumed  for  the  validity  of  statistical  theories  of  unimolecu- 
lar  reactions.  In  other  words,  the  probing  of  phase  space  by  the  dissociation  process  is  a 
coarse  grained  one. 

We  have  recently  argued30  that  more  than  the  completeness  of  the  sampling  of  phase 
space,  it  is  its  representativity  which  matters  to  insure  the  validity  of  statistical  theories.  It  is 
well  known46"58,59  that  the  individual  rate  constants  of  the  quasibound  states  in  an  energy 
interval  (£,  E  +  dE)  fluctuate  around  an  average  value  k(E )  which  can  be  estimated  by  the 
RRKM  theory.  Therefore,  by  representative  sampling,  we  mean  that  the  average  rate  calcu¬ 
lated  from  the  fraction  of  the  quasibound  space  which  has  been  sampled  agrees  with  the 
value  obtained  from  RRKM  theory.  In  the  case  where  the  localized  initial  excitation  sequen¬ 
tially  explores  regions  of  phase  space  coupled  to  different  product  channels,  selectivity 
might  be  expected.60 

We  have  studied30  the  competition  between  the  two  processes:  mode  mixing  and  disso¬ 
ciation.  Here,  the  two  relevant  observables  are  the  autocorrelation  function  IC(r)l2  which 
probes  the  persistence  of  the  memory  of  the  initial  excitation  and  the  survival  probability  of 
the  bound  subspace  B, 

P{t)  =  Tr[Bp(t)]  . 

The  entropy  5fl(r)  of  the  ME  density  matrix  normalized  in  bound  space  Bp(t)/P(t)  is  the 
measure  of  the  volume  sampled. 

Particular  attention  has  been  devoted  to  the  amount  of  mixing  of  the  initial  state  and  to 
the  role  of  overlapping  resonances.  The  fact  that  resonances  overlap  can  be  easily  under¬ 
stood  from  the  RRKM  equation  of  the  rate  constant:  the  ratio  of  the  width  fik  and  the  mean 
spacing  l/p(£)  is  large,  (ftp(E)k  =  N*(E  -  £0)»  1).  This  means  that  the  rate  matrix61  T  is 
off-diagonal  or,  in  classical  terms,  that  there  are  recrossings  at  the  barrier.  Or,  in  a  more  pic¬ 
torial  language,  there  is  a  healing62  of  the  broken  bond.  Under  such  circumstances,  the 
RRKM  rate  constant  is  only  an  upper  bound  and  the  system  might  have  enough  time  to 
sample  phase  space  before  dissociating.  Our  numerical  results30  show  that  due  to  recross¬ 
ings  the  rate  of  exploration  of  phase  space  is  indeed  increased  but  not  enough  to  lead  to  a 
representative  sampling  before  the  dissociative  process  occurs. 

Incomplete  mixing  and  sequential  exploration  of  phase  space  depend  in  a  very  sensi¬ 
tive  fashion  on  the  sharpness  of  the  localization  on  the  initial  state.  We  find  that  even  a 
small  amount  of  spreading  in  the  initial  state  (less  than  10%  of  the  available  phase  space  is 
initially  accessed)  is  sufficient  to  provide  a  substantial  sampling  of  phase  space  before  dis¬ 
sociation  and  thus  to  insure  the  validity  of  statistical  theories. 
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Appendix  C: 


INTRAMOLECULAR  DYNAMICS 

This  appendix,  by  F.  Remade  and  R.D.  Levine,  is  the  latest  review  of  our  work  on 
intramolecular  dynamics. 

It  was  presented  at  a  NATO  meeting  in  March  1992  and  is  published  as  a  chapter  in 
’Time  Dependent  Quantum  Molecular  Dynamics:  Experiment  and  Theory’. 
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INTRODUCTION 

The  temporal  evolution  of  an  initially  prepared  nonstationary  state  of  a  polyatomic 
molecule  is  discussed.  Particular  attention  is  given  to  the  short  time  evolution  which  can  be 
probed  by  spectroscopic  experiments.  It  is  suggested  that  the  time  evolution  is  often  sequen¬ 
tial,  with  separation  of  time  (or,  equivalently,  of  frequency  scales).  The  relevance  to  mode 
selective  chemistry  (in  both  the  gas  phase  and  in  solution)  is  emphasized.  In  particular  it  is 
argued  that  the  very  shortest  time  scale  identified  from  the  spectrum  is  typically  shorter  than 
the  rime  scale  for  mode  mixing  which  is  of  interest  in  statistical  theories  of  unimolecular 
dissociation. 

BACKGROUND 

On  the  one  extreme  are  direct  bimolecular  reactions1.  A  well  studied  prototype  is  the 
abstraction  of  an  H  atom  by  Fluorine2"4  RH+F  — >  R  +  HF.  Even  for  a  large  polyatomic, 
molecule  RH,  much  of  the  reaction  exoergicity  is  carried  away  as  internal  exdtation  of  HF. 
The  mechanical  picture  is  quite  clean  The  newly  formed  hot  HF  bond  is  coupled  to  R  for 
one  or,  at  most,  a  few  of  its  (very  rapid)  vibrations.  There  is  not  enough  time  for  the  energy 
in  the  ’hot  spot’  to  be  dissipated  over  all  the  available  degrees  of  freedom  of  RHF.  Our 
question  is  ’how  rapidly  does  such  a  hot  spot  delocalize’  and,  even  more  important  for  us, 
’is  the  dissipation  a  simple  process  or  does  it  have  fairly  well  defined  stages’.  That  there  is 
an  opposite  extreme  where,  given  enough  time,  the  initial  hot  spot  is  very  delocalized  is 
clearly  borne  out  by  the  chemical  activation  experiments5"10  especially,  those  of 
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Resonance  Raman  scattering  cross  sections,  as  a  function  of  the  excitation  fre¬ 
quency,  are  being  reported  for  excited  states  of  polyatomic  molecules  both  in  the  gas 
phase1-*  and  in  solution.5,6  The  time-dependent  interpretation7"9  of  Raman  spectroscopy 
provides  a  link  between  the  measured  frequency  dependent  cross  section  and  the  dynam¬ 
ics  on  the  excited  potential  energy  surface(s).  Specifically,  the  resonance  Raman  scatter¬ 
ing  amplitude  is  given  in  the  Condon  approximation  as  the  (half)  Fourier  transform  of  a 
time  dependent  amplitude  <  /I  exp(-iHt/H)\i  >  between  the  initial  non  stationary  state  on 
the  upper  surface  and  the  final  state.  This  enables  one  to  develop  both  physical  models 
and  semiquantitative  short  time  approximations  for  the  interpretation  of  the  main  features 
of  the  Raman  excitation  profile.  1,*-1S  As  discussed  in  detail  in  ref.  3  a  more  precise 
agreement  between  the  computed  time  crosscorrelation  function  and  the  observed  excita¬ 
tion  profile  requires  a  detailed  knowledge  of  the  relevant  excited  potential  energy  sur¬ 
face^). 

The  theoretical  route  up  to  now  has  been  from  the  potential  energy  surface(s)  via  a 
computed  time  crosscorrelation  function  to  the  spectrum.  The  half  (positive  times  only) 
Fourier  transform  relation  from  the  time  dependent  amplitude  to  the  scattering  amplitude 
precludes  the  direct  inversion  of  the  Raman  cross  section  to  yield  the  time  correlation 
function.  This,  in  contrast  to  absoiption  spectroscopy,  where  the  time  dependent  ampli¬ 
tude  can  be  directly  determined  as  an  inverse  Fourier  transform  of  the  frequency- 
dependent  absorption  cross  section.7,3,1^19 

The  ability  to  determine  the  time  autocorrelation  function  directly  from  the  observed 
absorption  cross  section  has  opened  the  way  towards  a  better  understanding  of  short  time 
intramolecular  dynamics.16-23  Even  when  the  potential  is  accurately  known  so  that  the 
time  autocorrelation  function  can  be  computed  (e.g.  ref.  19),  its  more  physically  intuitive 
meaning  suggests  that  it,  rather  than  the  cross  section,  provides  the  better  meeting  ground 
between  experiment  and  theory. 
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The  time  autocorrelatioa  function,  as  determined  from  the  absorption  cross  section 
gives  the  probability  to  remain  in  the  initially  accessed  non  stationary  state.  It  provides 
however  no  insight  as  to  where  the  initial  state  has  propagated  to.  The  latter  is  determined 
by  the  time  crosscorrelation  function.  The  inversion  of  excitation  frequency  dependent 
Raman  scattering  crc„s  sections  to  time  crosscorrelation  function  is  thus  of  general  inter¬ 
est  and  is  the  topic  of  this  Letter.  An  approach  which  is  also  capable  of  yielding  the  scat¬ 
tering  amplitude  but  which  is  quite  different  in  details  is  the  ’transform  theory’24  which 
seeks  to  determine  the  Raman  excitation  profile  from  the  absorption  cross  section. 

The  Raman  cross  section  for  scattering  from  state  i  to  state  /  is  proportional  to  the 
absolute  square  of  the  scattering  amplitude 

A 

afi(a))  =  (ifft)  jexp(icDt)Cfi(t)dt  .  (1) 

o 

Here  co  is  the  laser  frequency,  CfX[)  is  the  time  dependent  amplitude, 
<  /I  exp(-t//r/ft)lf  >,  we  are  after  and,  for  the  time  scale  of  interest  to  us  here,  we  did  not 
include  in  (1)  the  radiauve  damping  of  the  excited  state.  There  are  two  technical  difficul¬ 
ties  in  inverting  (1).  The  first  is  that  the  experiment  only  determines  la/-,!2,  whereas  for 
inverting  (1)  one  needs  the  real  and  imaginary  parts  separately.  The  half  Fourier  trans¬ 
form  in  (1)  suggests  that  be  re6.  led  as  the  full  Fourier  transform  of  the  ’causal’ 

function  U{t)Cfi{t)  where  U(t)  is  the  unit  step  function.  Most  texts  which  discuss  the 
Fourier  transform  of  causal  functions  explicitly  or  implicitly  make  the  assumption  that 
they  are  real  and  proceed  to  establish  a  Hilbert  transform  relation  between  the  real  and 
imaginary  parts  of  their  Fourier  transform,  here  or^fry).  However,  our  C} -,(t)  is  not  neces¬ 
sarily  real  so  that  all  that  one  can  establish  in  general  is25  that  the  even  and  odd  parts  of 
afi(a))  are  related  as  Hilbert  transforms.  We  therefore  proceed  as  in  the  derivation  of  dis¬ 
persion  relations  in  other  branches  of  scattering  theory26  and  assume  that,  as  a  function  of 
s  =  id),  (Zfi(s)  is  analytic  and  has  no  zeroes  in  the  right  half  of  the  complex  s  plane, 
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Rc  J>0.r  This  implies  that  ln(a^(^))  is  analytic  for  Re  s  >0  which  is  the  condition  we 
need  for  a  unique  inversion.  One  can  also  use  this  condition  to  establish  dispersion  rela¬ 
tions  between  the  real  and  imaginary  parts  of  ln(ccji(s)).  Since  ReIn(<Zy;)  =  ln(!ay;l)  is 
known  from  experiment  these  dispersion  relations  can  in  principle  be  used  to  determine 
the  phase  of  a ft  (which  is  the  imaginary  part  of  ln(ar^)).  We  ha\e  implemented  such  a 
procedure  and  will  report  on  it  in  detail  in  the  full  account  of  this  work.  Here  we  just  say 
that  this  route  is  workable  and  quantitative  but  not  quite  as  numerically  accurate  as  the 
route  we  now  describe. 

Specifically,  we  write 

la ys(<»)l2  =  1/  2  A,  exp (io)tr)  .  /•  (2) 

r»-« 

The  2n  + 1  expansion  coefficients  Ar  satisfy  A_,  =  A*  because  lor/2  is  real.  The  time 
points  r,  are  not  necessarily  all  equally  spaced.  The  numerical  values  of  both  the  A/s  and 
the  r/s  is  determined  by  a  procedure  of  maximum  entropy28,29  as  is  discussed  below.  The 
assumed  analytical  properties  of  afi(o))  are  needed  to  show  (e.g.,  A.  Papoulis,  ref.  27, 
loc.  cit .)  that  (2)  leads  to  a  unique  form  a^a)  which  we  write  as 

n 

afi{a))  =  l/(£(ar  +  z&,)  exp(/ryr,))  with  real  a’s  and  b’s.  The  A/s  are,  of  course,  bilinear 

r-0 

in  the  a/s  and  br’s.  The  condition  A*  =  A_r  insures  that  the  number  of  unknowns  be  the 
same.  IF  we  take  all  the  r/s  to  be  equally  spaced  then  (2)  is  equivalent  to  the,  so  called. 
Burg  spectrum30,3 1  and  the  a’s  and  b’s  can  then  be  determined  by  the  Levinson  algo¬ 
rithm.31  This  provides  an  easy  route  for  an  initial  guess  and  does  result  in  a  realistic 
Cfi(t).  However,  it  leads  to  a  periodic  function  la^(o;)l2,  which  is  not  physical,  and  a 
superior  representation  of  afi((o)  is  obtained  when  the  r/s  are  not  all  equally  spaced. 

From  the  experimental  spectrum  one  can  determine  the  Fourier  transform  cf 
The  functional  form  (2)  is  obtained2*  as  that  of  the  most  probable  (for  a  distri¬ 
bution  of  maximum  entropy)  form  of  lar^(ry)!2  subject  to  given  numerical  values  of  /?^(r) 
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at  the  given  set  of  times  f,.  A0  is  the  multiplier  which  insures  that  the  area  under  lafi(co)t2 
equal  Rfi( 0).  The  values  of  the  Lagrange  multipliers,  Xr,  are  determined  by  the  algorithm 

of  ref.  32  with  two  modifications:  (1)  The  entropy  is  written  as  J  dco  ln!or^(<y)l2  with  the 

maximum  being  sought  not  directly  as  a  function  of  the  A,’s  but  by  varying  the  (equiva¬ 
lent)  set  of  unknown  ar' s  and  br's  on  which  itself  depends.  (2)  Once  a  given  value 
of  the  maximum  of  the  entropy  is  secured,  by  optimising  the  ar’s  and  b,'s  for  a  given  set 
of  r,’s,  a  new  set  of  tr's  is  sought  such  that  the  new  maximal  value  of  the  entropy  is 
lower.  Elsewhere,  we  shall  describe  this  min-max  aspect  of  the  algorithm  in  detail.  In 
general,  we  start  with  equally  spaced  tr' s  and  and  with  trial  values  for  the  ar's  and  br's 
determined  using  the  Levinson  algorithm.31  The  time  dependent  amplitude  Cfi(r)  is 
obtained  as  a  Fourier  transform  of  afi{co).  Special  care  is  needed  to  obtain  the  overall 
phase  of  which  depends  on  its  asymptotic  behavior.33  Even  when  this  care  is  exer¬ 

cised,  a  use  of  a  library  fast  Fourier  transform  algorithm  will  often  spoil  it.  A  quick  but 
risky  route  is  to  take  advantage  of  C ^(f)  being  guaranteed  causal  and  shift  the  zero  of 
time  to  where  C fi(t )  starts  to  differ  from  zero. 

Figure  1  compares  three  time  crosscorrelation  functions  \Cfi(t)l2  computed  for 
iodobenzene  using  equations  (3)-(6)  of  ref.  4  with  those  obtained  by  the  present  inversion 
from  \afi{co)\2.  The  three  results  shown  correspond  to  a  final  state  with  1,2  and  3  vibra¬ 
tional  quanta  in  the  vu  mode  of  the  B  electronic  state.  Note  how  the  method  keeps  track 
of  the  systematic  changes  in  the  time  crosscorrelation  function  for  different  final  states. 
Nine  time  values  (n  =  9  in  equation  (2))  were  used  to  obtain  the  results.  We  do  not  show 
a  comparison  of  the  correct  excitation  profile  and  equation  (2)  because,  with  n  =  9  and 
the  min-max  algorithm,  the  fit  is  quite  tight. 

General  properties  of  the  Fourier  transform  determine  the  requirements  on  the 
experimental  data  for  a  successful  inversion.  In  particular,  one  needs  a  spectral  range 
broad  enough  to  generate  the  entire  envelope  t  f  the  time  crosscorrelation  function.  Also, 
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a  good  reproduction  of  the  very  short  time  dynamics  requires  an  accurate  representation 
of  the  asymptotic  (high  and  low  oS)  shape  of  afi(co).  This  is  noticeable  in  the  bottom 
panel  of  the  figure  for  the  0  -»  1  scattering  cross  section  where  the  asymptotic  1  fco 
dependence  of  is  not  well  reproduced  since  we  intentionally  truncate  our  input  lor^I2 
to  mimic  realistic  data.  On  the  other  hand,  note  the  good  reproduction  of  the  relative 
heights  and  positions  of  the  main  features. 

Our  considerations  can  also  be  used  for  the  absorption  cross  section.  Equation  (2) 
determines  both  the  position  and  the  intensities  of  all  transitions.  It  is  thus  of  considerable 
interest  in  connection  with  the  distribution  of  spectral  spacings  and  intensities.  Elsewhere 
we  shall  show  how  the  present  developments  can  be  viewed  as  a  direct  generalization  of 
the  maximum  entropy  approach34  to  these  issues. 


It  has  been  shown  possible  to  invert  frequency-dependent  Raman  scattering  cross 
sections  so  as  to  determine  time  crosscorrelation  functions.  A  new  maximum  entropy 
algorithm  provides  an  essentially  quantitative  inversion  and  an  accurate  fit  of  the  cross 
section  itself.  Only  the  bare  essentials  of  the  approach  have  been  described  here  and  a 
full  account  is  in  preparation  for  publication.  Other  aspects  not  discussed  in  detail  include 
the  implications  for  the  analysis  of  the  exploration  of  phase  space  and  the  representation 
of  spectral  fluctuations. 
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Figure  legend 

Tlie  directly  computed  time  crosscorrelation  functions  IC^r)!2  (left  column)  and  the 

results  of  the  inversion  using  \afXco)\z,  (right  column)  vs.  time  (in  fs).  Other  details  in  the 
text. 
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